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^ | . Abstract 

■ The present paper gives an overview of the recent developments in the descrip- 

tion of critical behavior for Hamiltonian perturbations of hyperbolic and elliptic 
systems of partial differential equations. It was conjectured that this behavior can 
, be described in terms of distinguished Painleve transcendents, which are univer- 

sal in the sense that they are, to some extent, independent of the equation and 
<-j , the initial data. We will consider several examples of well-known integrable equa- 

tions that are expected to show this type of Painleve behavior near critical points. 
The Riemann-Hilbert method is a useful tool to obtain rigorous results for such 
equations. We will explain the main lines of this method and we will discuss the 
universality conjecture from a Riemann-Hilbert point of view. 



1 Introduction 

m 

IT) 

■ We will discuss the asymptotic behavior of solutions to Hamiltonian perturbations 

of hyperbolic and elliptic systems of partial differential equations near critical points 
of the unperturbed system. We will give an overview of recent developments in this 
area using the Riemann-Hilbert (RH) approach. The systematic study of the critical 
behavior of such Hamiltonian systems was initiated by Dubrovin and collaborators 
[TT1 [J5J [T9l [201 EE] i n a series of papers where it was conjectured that the local behavior 
^ ■ of solutions near critical points is to a large extent universal, and that it can be described 

in terms of distinguished Painleve transcendents. The conjectured behavior has been 
proved and verified numerically in various special cases. We will focus mainly on 
particular cases of Hamiltonian PDEs or systems of PDEs that can be written in a 
so-called Lax form and that consequently can be solved using the direct and inverse 
scattering transform, which can be formulated as a RH problem. The Deift/Zhou 
steepest descent method then provides the tools to obtain asymptotics for the RH 
problem and for the solution to the system of equations. 

Let us first consider the simple quasi-linear hyperbolic equation ut + 6uu x = 0, 
which is known as the Hopf equation or in-viscid Burgers' equation. It is a classical 
fact that solutions to this equation exist only for small times t > and develop shocks 
at a certain time. Given initial data uq(x), the solution is given by the method of 
characteristics in the implicit form 

u(x,t) = u (0, x = 6tu (0+C- (1-1) 
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Figure 1: The solution to the Hopf equation with initial data uo(x) = — sech 2 (x) for 
t = 0, t = 0.11, and t = 0.216506. 



For simplicity, we will only consider smooth initial data that decay sufficiently rapidly 
at ±oo. One observes that the derivative of the solution blows up at the critical time 

1 

max£ eK [-6u' (£)]' 

see Figure [TJ We write x c for the point where the derivative of u(.,t c ) blows up, and 
u c — w(x c ,t c ). 

The blow-up of the first derivative is called the gradient catastrophe, and we say 
that x c and t c are the point and time of gradient catastrophe. Locally near x c , the 
solution at the critical time behaves (for generic initial data) like 

u(x,t c ) = u c - c(x - x c ) 1//3 + 0(x - x c )^ 3 , as x — > x c . (1.2) 

In order to avoid the gradient catastrophe, one can perturb the Hopf equation. The 
Burgers' equation u t +6uu x — e 2 u xx = is an example of a dissipative regularization, but 
we will concentrate on dispersive or Hamiltonian perturbations of the Hopf equation. 
In [17] , all Hamiltonian perturbations to the Hopf equation ut + Quu x = up to order 
4 have been classified, which lead to the class of equations 

+2p' (5u xx u xxx + 3u x u xxxx ) + p"{7u x u 2 xx + §u 2 u xxx ) + 2p"' u x u xx = 0, (1.3) 

where p = p(u) and c = c{u) are arbitrary smooth functions. For c = 12, p = 0, we 
have the KdV equation 

u t + 6uu x + e 2 u xxx = 0. (1.4) 

For sufficiently smooth initial data uq(x) = u(x,0,e), solutions to the KdV equation 
exist for all times t > 0. For small times, u and its derivatives are bounded and the 
third derivative term will only contribute up to order C(e 2 ) for small e. Near the time 
and point of gradient catastrophe for the Hopf equation, solutions to KdV start forming 
oscillations, and the x-derivative of the solution is no longer uniformly bounded. The 
contribution of the dispersive term e 2 u xxx in the equation becomes visible in this region: 
for t slightly bigger than t c , it is well-known that there is an interval [x~ (t) , x + (t)} 
where the KdV solution u(x, t, e) shows oscillatory behavior. Outside this interval, the 
KdV solution can still be approximated by a continuation of the Hopf solution. The 
oscillations have been studied extensively in many works and are given asymptotically 
in the small dispersion limit in terms of the solutions to the Whitham equations for KdV 
and in terms of Jacobi elliptic functions and elliptic integrals [121 113 E21 [391 El EH] ■ 
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At least before the gradient catastrophe, the small e behavior for general equations 
in the family (jl.3p is expected to be similar to the behavior for KdV, but few analytical 
results are available. Even existence of solutions with smooth initial data for small 
times has not been established in general. For small times, it is not hard to believe 
that the contributions of the terms of order e 2 and e 4 will be small. One expects that 
u(x,t,e) = u(x,t) + C(e 2 ) for t < t c — 5, where u(x,t) is the solution to the Hopf 
equation with initial data uq(x). It is less obvious what happens for x and t near the 
point and time of gradient catastrophe. This problem was addressed by Dubrovin in 
|17j . where universality of critical behavior for Hamiltonian perturbations of hyperbolic 
PDEs was conjectured. 

Universality conjecture: the hyperbolic case 

It is expected that a generic solution to an equation of the form (jl.3p has an expansion 
of the form 

u(x, t, e) = u c + de 2 / 7 Ef( CB a ~ gc ~; 3(t ~' c) , c 4 ^) + 0{e^), (1.5) 

as e — > and at the same time x — > x c and t — >• t c at appropriate speeds, where x c and 
t c are the point and time of gradient catastrophe. The constants 01,02,03, C4 depend on 
the equation and the initial data, but not on x, t. The function U(X, T) is expected to 
be universal, i.e. independent of the equation and independent of the choice of initial 
data. It should be a solution to the fourth order ODE 

X = TU- Qf/ 3 + ±(U X + 2UU XX ) + ^Uxxxx^j (1.6) 

which is smooth and real for all real values of X and T. The existence and uniqueness 
of such a solution was part of the conjecture, but the existence of a real smooth solution 
with the asymptotic behavior 

U{X,T) = t{§\X\) 1,z + 0{\X\- 2 ' Z ), asl^oo, (1.7) 

for any fixed TsK, has been proved in [8]. The smoothness of the Painleve transcen- 
dent U{X,T) was already conjectured long time ago for T = 0H2], and appeared 
later to be related to the Gurevich-Pitaevskii special solution to KdV [321 113 US] . In 
the physics literature, U(X, T) shows up in the study of ideal incompressible liquids 
[371 1381 [29] and in quantum gravity [16 1 . It appears also in the description of the lo- 
cal eigenvalue behavior in unitary random matrix ensembles with singular edge points 
019]. 

Substituting the asymptotics (|1.7p in (|1.5p for t = t c , we recover (jl.2p in the limit 
e — > 0. Equation (]1.6p is known as the second member of the Painleve I hierarchy and 
has, given T G C, solutions that are meromorphic in X with an infinite number of 
complex poles. The smoothness of U means in other words that no poles lie on the real 
line. It is remarkable that U(X,T) is an exact solution to the KdV equation written 
in the form U T + UU X + j^Uxxx = 0. 

Loosely speaking, the conjecture suggests that the local behavior of solutions to 
Hamiltonian perturbations of the Hopf equation is universally described in terms of 
U(X,T): the same Painleve transcendent appears independent of the initial data and 
independent of the equation. U can be seen as the function which describes the transi- 
tion between the region where the small dispersion asymptotics are determined by the 
Hopf equation and the region of oscillatory behavior. 
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Critical behavior in terms of U(X, T) is expected [18J for more general Hamiltonian 
perturbations of strictly hyperbolic systems of equations of the form 



u t + A(u)u x = 0, 



u 



(u{\ 
\u n ) 



(1.8) 



where A(u) is a n x n matrix-valued function of u, strict hyperbolicity means that A 
has real and distinct eigenvalues in the domain of the (x, t)-plane under consideration. 

In the case of the KdV equation and KdV hierarchy (see below), the universality 
conjecture has been proved for a class of negative, real analytic initial data u${x) with 
one local minimum, and which tend to rapidly as x — > ±oo. For the KdV equation, 
we have fjl .5j) with the values [6] 
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ci 



(8fc) 2 / 7 ' 



C-2 



(8k) 1 / 7 ' 



c 3 = 6u c , c 4 



12 



(8k) 3 / 7 ' 



where k = — f'£'(u c ), with /l the inverse of the decreasing part of the initial data uo. 
The expansion holds in a neighborhood of (x c , t c ) which shrinks with e, to be more 
precise it holds in the double scaling limit where e — > 0, x — > x c , t — > t c in such a way 
that 



lim C2 



%c C 3 (t t c 



6 

er 



X, lim C4 - 



t -t c 

A/7 



T, 



(1.9) 



Examples 

We will now list a number of equations and systems of equations that are expected to 
belong to the universality class for which U(X, T) describes the solutions locally near 
the gradient catastrophe. 

(i) The m-th time flow of the KdV hierarchy is given by 

Ut - (-i) m d x c m = o, (i.io) 

where C m is the Lenard-Magri recursion operator defined by 

d x C m = (e 2 <9^ + Aud x + 2u x ) £ m -i, C = u. (1.11) 

For m = 1, we have the KdV equation, and if we impose C m to be zero for u = 0, 
the second and third equation in the hierarchy are given by 

ut — 2>Qu 2 u x — e 2 (20u x u xx + Wuu xxx ) — e 4 u xxxxx = 0, (1-12) 

u t + U0u 3 u x + e 2 (70u 3 x + 280uu x u xx + 70u 2 u 3x ) 

+ e 4 (70u 2x u 3x + A2u x U4x + 14uu 5x ) + e 6 u 7x = 0, 

(1.13) 

where we have written Uj x for the j-th partial derivative of u with respect to 
x. The universality conjecture has been confirmed numerically for the second 
member of the hierarchy in [22J, and an expansion of the form (11.50 was proved 
analytically in [7] for all equations in the hierarchy. 
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(ii) For the Camassa-Holm (CH) equation 

u t +(3u+2n)u x -e 2 (u xxt + 2u x u xx + uu xxx ) = 0, t > 0, x £ R, k G R, (1.14) 

numerical results have been obtained in [30] supporting the universality conjec- 
ture, although the equation is not precisely of the form (jl.3j) or (11.8P , see [31] . 

(iii) The de-focusing nonlinear Schrodinger (NLS) equation 

e 2 

ie^ t + -^ xx - |Vf V = 0, (1.15) 



can be transformed to the system 

Ut + vu x + uv x = 
v t + vv x -Ux + ^(jz^? - 



(1.16) 



where u = — \ip\ 2 , v = ^ — ^f-J. This is a Hamiltonian perturbation of the 
system 

:), + (-. :)©.-* (1 - 17) 

which is hyperbolic because the eigenvalues of the coefficient matrix are v db \f—u. 
(iv) For the Kawahara equation 

ii t + 6im x + e 2 u xxx = e 4 -u xx:r:rx (1-18) 
and generalized KdV equations 

u t + 6u n u x + e 2 u xxx = 0, (1-19) 
the universality conjecture has been verified numerically in |22] , 

Universality conjecture: the elliptic case 

The focusing nonlinear Schrodinger equation 

ieifjt + -i) xx + = (1.20) 

can be seen as a Hamiltonian perturbation of an elliptic system, since it is transformed 
to the system (|1.16p after the substitutions u = \ijj\ 2 : , v = ^ ( % — This is a 



2i \ 1/1 i/j* 

Hamiltonian perturbation of the system (|1.17p which is then of elliptic type since the 
eigenvalues v ± iy/u are complex conjugate for positive u. The Cauchy problem for 
this equation is ill-posed and the blow-up phenomenon is essentially different than for 
hyperbolic systems. The system (|1.17p can be solved using a hodograph transform 
method up to a critical time t c where the spatial derivative of u blows up at a point 
x c . This happens when the peak at a local maximum of the solution becomes more 
and more narrow (or focusses). This critical point is called a point of elliptic umbilic 
catastrophe, and a conjecture has been formulated in [21j concerning the local behavior 
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of the NLS equation near this point. The key role in this conjecture is played by the 
Boutroux tritronquee solution to the first Painleve equation 

Qzz = Z + 6Q 2 , (1.21) 

characterized by the asymptotic behavior 

f~Z 

Q(Z) J-, as Z -> ±00 in the sector | argZ| < ^. (1.22) 

It is known that this is a meromorphic function with infinitely many poles in the sector 
of the complex plane where | argZ| > In the complementary sector | axgZ\ < 
Q has no poles for sufficiently large Z (this would contradict the asymptotic behavior 
(I1.22D ) and no poles at all on the positive half- line [33], but this does not exclude the 
possibility that a finite number of poles is present in this sector. It was conjectured by 
Dubrovin, Grava, and Klein in |21j that Q has no poles at all in the sector | arg Z\ < 
The second part of the conjecture, in analogy to the hyperbolic case, describes 
the behavior of the solution to the focusing NLS equation: it is expected to have 
an asymptotic expansion where the first error term to the unperturbed solution is of 
order e 2 / 5 and can be described completely in terms of the tritronquee solution Q(Z) 
evaluated at complex arguments Z. 

Substantial progress has been made on a proof of the second part of the conjecture 
in the recent work [2], we will comment in more detail on this later on. The poles of 
Q have been studied in [3D] , but the first part of the conjecture is still open. It was 
discussed in [18] that an expansion in terms of Q(Z) is not exclusively expected for the 
focusing NLS equation, but also for more general Hamiltonian perturbations of elliptic 
systems. The tritronquee solution Q is also relevant for the asymptotics of recurrence 
coefficients of certain orthogonal polynomials on a complex contour [23\ I25j. 

Outline 

We will not discuss classification results or try to formulate the most general form 
of the conjectures, for these topics we refer to [T71 [181 032 023 (ZD E2]- Rather we 
want to focus on techniques that can be used to prove the universality conjectures in 
various special cases. The best one can do so far is to try to prove the universality 
conjectures for systems of equations that can be solved using the inverse scattering 
transform. Then one can formulate a RH problem which characterizes solutions to the 
system of equations, and an asymptotic analysis of this RH problem can in principle 
lead to the critical asymptotics in terms of the Painleve transcendents U and Q, which 
themselves can be characterized in terms of a so-called model RH problem. Without 
going into the technical details, we will give an overview of this method in the case 
of the KdV equation (following [12\ \T3[ [6] ) , and indicate the differences for the KdV 
hierarchy (following [7]). Afterwards we will recall the RH problems associated to the 
Camassa-Holm equation [4| and the de- focusing and focusing NLS equation [44], and 
explain the main lines of the procedure that can be followed to obtain asymptotics 
for solutions of those equations. We will also give heuristic arguments supporting the 
universal nature of the critical behavior from the Riemann-Hilbert point of view. 
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2 RH problems for Painleve transcendents 

The Painleve transcendents U and Q, which are the central objects in the universality 
conjectures, can be characterized in terms of RH problems. We will state the RH prob- 
lems here, and we will explain that the absence of poles is equivalent to the solvability 
of the RH problems. 

2.1 RH problem for Pf solution U 

The function U (X, T) is characterized by the following RH problem. 



RH problem for ^: 

(a) # = V((;X,T) is analytic for ( G C\T, with T = lUe^l + Ue"TM+ oriented 
as in Figure [21 

(b) ^ has continuous boundary values ^± (+ referring to the left side of the contour, 
— to the right side according to the orientation in Figure [2]) on T \ {0}, and we 
have 

* + (C) = *-(C)(_° 1 J), forCe IT, (2.1) 

*+(C) = *-(C)(j J), forCGM+, (2.2) 

* + (C) = *-(C)(j J), forCGe ± ¥R+. (2.3) 

(c) has the following behavior at infinity, 

*( C ) = fl+ MX,T) + 0{{ --2^ c ^ Ne -0^,T)^ t (2 .4) 

where 03 = ( J J , and N and 8 are given by 



N 



^(_\ ^(C;^T) = I ^C 7/2 -|c 3/2 + ^C 1/2 , (2.5) 



and where A\(X,T) is a matrix independent of £• 

It was proved in [8] that this RH problem has a solution for all real values of X, T, and 
that the function 

U(X,T) = 2A 1:11 (X,T) - A\ l2 {X,T) (2.6) 

is a real pole- free solution to equation (|1.6p which has the asymptotic behavior (|1.7j ). 
General solutions to (|1.6|) are characterized by a similar RH problem, but with jumps 
also on the four rays e^^TR -1- and e^^nR -1- [33]. The jump matrices are all triangular 
with ones on the diagonal, and the triangular structure is as follows: the jump matrices 
are upper-triangular on R, e 7 ~R + , and lower-triangular on e 7 M and e 7 K . 
The off-diagonal entries of the jump matrices are called Stokes multipliers and satisfy 
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Figure 2: The jump contour T and the jump matrices for 



certain constraints. Given T E C, there is a one-to-one correspondence between the 
set of admissible Stokes multipliers (also called the monodromy surface) and the solu- 
tions of the Pf equation. The RH problem is known to be solvable if and only if the 
corresponding Painleve solution is analytic at the value (X, T), or in other words the 
RH problem is not solvable if and only if (X, T) is a pole of the corresponding Painleve 
solution. 

2.2 RH problem for PI 

The Boutroux tritronquee solution Q(Z) of the Painleve I equation is characterized by 
another RH problem. 

RH problem for $: 

(a) <£> = Z) is analytic for (£C \ V, with f = tU e^R + U e" 2 f i R+ oriented 
as in Figure [2721 

(b) 3> satisfies the following jump relations on T' \ {0}, 

M0 = M0(- o)' forCeM", (2.7) 

MO = Mo(- J), forCGM+ (2.8) 

MO = MO (J f), forCGe ±i "M+. (2.9) 

(c) $ has the following behavior at infinity, 

m = {i+^f 1 + o(c 2 ))c^^(\ \y-^ z ^, (2.io) 

where a is given by 

a{C-Z) = ^/ 2 -ZC 1 ' 2 . (2.11) 
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Figure 3: The jump contour V and the jump matrices for <3? 



The function 

Q(Z) = Bf 21 - 2B h22 (2.12) 

is the tritronquee solution to Painleve I [35] , and given Z £ C, the RH problem for 
<3? is solvable if and only if Z is not a pole of Q. A general solution to the Painleve 
I equation can be obtained in terms of a RH problem with triangular jumps on the 
contour KUe ±! ^M + Ue ±i TK + [351 123] . 

In order to discuss the solvability of this RH problem, let us first introduce the 
homogeneous RH problem for <3?: we say that <3?o is a solution to the homogeneous RH 
problem if it satisfies conditions (a) and (b) of the RH problem for <3>, and it has the 
asymptotic behavior 

= o(£-3/4) e -a(C;Z)<T 3) asC^oo. (2.13) 

Following the vanishing lemma procedure [24[ [26] [27] , the solvability of the RH problem 
for $ is equivalent to the fact that the homogeneous version of the RH problem has 
only the trivial zero solution: the RH problem for <3? with parameter Z E C is solvable 
if and only if the homogeneous RH problem has only the vanishing solution <I>o = 0. To 
prove this, one can follow the procedure developed in |10[ [26j [27] (which works for a 
fairly general class of RH problems) and show that a certain singular integral operator 
related to the RH problem is a Fredholm operator of index zero. 

So in order to have smoothness of Q for | argZ| < it suffices to show that the 
homogeneous RH problem has only the zero solution (such a result is called a vanishing 
lemma) for | argZ| < 4p Then the RH problem for <I> is solvable, and Q has no pole 
at Z. 

A vanishing lemma has been proved for the Pf RH problem for real X and T jS], 
but this turns out to be considerably more complicated in the PI case, even for positive 
real values of Z. 

3 Small dispersion asymptotics for solutions to KdV: the 
Riemann-Hilbert approach 

We will now give an overview of the RH method to obtain asymptotics for solutions to 
the KdV equation in the small dispersion limit. The method consists of three steps: first 
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one has to construct a RH problem that characterizes solutions to the KdV equation. 
This can be done using scattering theory for the Schrodinger operator. The jump 
matrices in the RH problem will depend on the initial data via a reflection coefficient. 
The second step requires asymptotic information about the reflection coefficient as 
e — > 0. Such results were obtained in [33] using WKB techniques. Once this information 
is gathered, one can start with the final step which is the asymptotic analysis of the 
RH problem using the Deift/Zhou steepest descent method 



3.1 Construction of a RH problem related to the Schrodinger opera- 
tor 

The key observation here is that the KdV equation can be written in the Lax form 
Lt = [L,A], where L is the Schrodinger operator, and A is an antisymmetric operator: 

Lf = e 2 8 2 x f + uf, Af = A = 4e 2 d 3 J + 3(ud x f + d x (uf)). (3.1) 

We consider potentials u{x) that are smooth and decay sufficiently fast at ±oo. For 
a general real potential u(x), L has a finite number of purely imaginary eigenvalues 
Xj for which the linear Schrodinger equation Lf = Xf has solutions decaying both at 
+oo and at — oo. For A < 0, it admits solutions ip± and (f>± satisfying the oscillatory 
asymptotic conditions 

lim V±(A;x,i,e)e ± ^ v/rix = 1, (3.2) 

X— >+oo 

lim 4> ± (X;x,t,e) Tl ^ x = 1. (3.3) 

x— >— oo 

Those solutions are called Jost solutions. They are related by a connection matrix 
given by 

*-W> = $A;M)) ' (3 - 4) 

The ratio r(A;i,e) = -(X;t, e) is called the reflection coefficient for the Schrodinger 
equation. 

It is well-known that ip+ (?/>_) and 0+ (0-) can be extended analytically to the 
upper (lower) half plane. The vector-valued function 



M(X;x,t,e) 



(3.5) 

(^V- 0-) e^ ( - A)1/2<T3 , asAeCT, 



is meromorphic in C \ R and it can be showed using (|3.4|) and asymptotic properties of 
the Jost solutions that M satisfies the following RH conditions [441 fTTj [T2l [6] 

RH problem for M 

(a) M(A; x, t, e) is meromorphic for A G C\K, 
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(b) M has continuous boundary values M±(X) for A6l\ {0} that satisfy the jump 
conditions 



M+(A) = Af_(A) 



1 r(A;t,e)e^ :E (- A ) 1 
^(A;t,e)e^-(-^ 1/2 1 - |r(A; t, e)| 2 



M+(A) = M_(A)<7i, cji 



1 

1 



for A < 0, 

(3.6) 

for A > 0. 

(3.7) 



(c) As A — > oo, we have M(X; x, t, e) — > (l l) . 
The potential u(x,t, e) can then be recovered from the equation 

u(x,t,e) = -2ie lim ^T^Xd^Mu^x^t^e) - 1). (3.8) 

A— >oo 

In general M has a finite number of simple poles on the imaginary axis, which are the 
eigenvalues of the Schrodinger operator. One has to impose additional conditions on 
the residues of M to have a unique RH solution. If the potential u{x) is negative, the 
Schrodinger equation has no point spectrum and M is holomorphic in C \ K. Then one 
can also use the x-derivatives of the Jost solutions to construct the second row of a 
2x2 matrix M satisfying the same jump condition, and the asymptotic condition 

M(A) ~ ( . I— . l r — ) , as A ^ oo. (3.9) 



-A 

A matrix RH problem is for several reasons more convenient to work with than a vector 
RH problem. 

The above construction works for any potential u(x) decaying sufficiently fast at 
±oo, it is not necessary that u = u(x,t,e) solves the KdV equation. However, for 
general potentials u(x,t, e), the time-dependence of r is not explicitly given. On the 
other hand if u(x,t,e) satisfies a suitable Lax type equation Lt = [L,A], the time 
evolution of r can be explicitly calculated, and the eigenvalues Aj are independent of t. 
In the KdV case, one has [28] 

r(A;t,e)=ro(A;e)eT*(-^) 1 . (3.10) 

This is a crucial observation: if we know the initial reflection coefficient corresponding 
to the initial data uq(x) (the direct scattering transform), then we know the reflection 
coefficient at any time t > 0. In order to find the solution u(x, t, e) at time t > 0, it 
suffices to solve the RH problem for M at time t > (the inverse scattering transform) . 



3.2 Semi-classical asymptotics for the reflection coefficient 

The solution to the Cauchy problem for KdV is characterized in terms of the RH 
problem for M. The jump matrix for M depends on the reflection coefficient ro, which 
depends in a simple explicit way on t, but in a much more complicated way on the 
spectral variable A and on e. Since we are interested in the small e behavior of the RH 
solution M, inevitably we need small e asymptotics for the initial reflection coefficient 
tq. Detailed asymptotics were obtained in |43| for negative initial data with a single 
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negative hump that are analytic in a sufficiently large neighborhood of the real line, and 
we will concentrate on such data in what follows. For A < —1, ro(A;e) decays rapidly 
as e — > 0, and for —1 < A < 0, it is oscillatory and has the leading order behavior 

ro(A; e) ~ i exp ^— — p(X)\ , as e — > 0, (3-11) 

where p(X) = \ ^t|~~x^£> anc ^ II is the inverse of the decreasing part of the initial data 
no- Near —1, there is a subtle transition between the decay and oscillatory behavior. 

Substituting tq in the jump matrix of the RH problem by its leading order asymp- 
totics (13. llh . we obtain the jump relation M + (A) = M_(X)vm(X) for A 6 M, with 

vm(X) ~ I for A < —1, 

(3.12) 

m ( 1 «exp(-f (p(A)-a(A)))\ , 

w "W~^a«p(S(p(A)-a(A))) 1 <V W J for-KA<0, 

(3.13) 

vm(X) = o\ for A > 0, 

(3.14) 

where a(X) = x(— A) 1 / 2 + 4i(— A) 3 / 2 . It is not obvious at all that one can proceed 
with the RH analysis, neglecting the corrections to the leading order behavior of the 
reflection coefficient. This needs to be justified, and especially for A near —1, technical 
issues need to be handled. We will not address this here, since it is our goal to give a 
flavour of the method only, we refer the reader to [6j [7] for details. 

3.3 Steepest descent analysis of the RH problem 

The Deift/Zhou steepest descent method was developed in [14] and applied to the 
small dispersion limit for the KdV equation in [12\ 113] . The method consists of a 
number of invertible transformations of the RH problem, and the goal of this series of 
transformations is to obtain a RH problem for which the jump matrices are uniformly 
close to the identity matrix as e — > (or in a more complicated double scaling limit 
where x — > x c , t — > t c simultaneously with e — > 0), and with a solution that tends to 
/ as A — > oo. For such RH problems, one can apply small norm theory to prove that 
the solution is uniformly close to the identity matrix for small e and to obtain the next 
terms in the asymptotic series of the solution. Inverting each of the transformations, one 
obtains asymptotics for M from the asymptotics for the final small norm RH solution. 

The first transformation requires the construction of a (^-function, which has suitable 
jump conditions and asymptotic behavior. For x,t sufficiently close to x c ,t c , the g- 
function g = g(X; x,t) is defined by the scalar RH conditions 

(a) g is analytic in C \ [u, +oo), 

(b) g has the jump relations 

5+ (A) + 5 _(A) =0, for A G (0, +oo), (3.15) 

g+(\) +g-(X) - 2p(\) + 2a(A) = 0, forAe(u,0), (3.16) 

(c) g(X) = ©(A" 1 / 2 ) as A -> oo. 
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Given it, this scalar RH problem has a solution (which is unique if one imposes suitable 
behavior of g near u and 0) that can be expressed explicitly as an integral transform of 
the inverse of the initial data. In order to analyze the KdV solution in the vicinity of 
x c and t c , the choice of the point u is crucial: it has to be chosen as u = u c = u(x c , t c ). 

The g-function is an essential ingredient for the asymptotic analysis of the RH 
problem for M. It clears the road for a transformation of the RH problem to one with 
jump conditions of the following form, on a jump contour as the one shown in Figure 

m 



RH problem for S 

(a) S is analytic in C \ £5, 

(b) S+(A) = S-(X)vs for A G Eg, and as e — > 0, the jump matrix vs has the leading 
order behavior 



on Ei, 



v s (X) ~ < 



'(: 






1 




1 






(- 






i 







on E 2 , 

as A G (n c , 0), 
as A G (-00, -1 - 61) U (0, +00). 



(c) 5(A) = (J + ©(A- 1 )) M-^ _ as A -> 00. 

The function <j> appearing in the jump matrices is given by 

4>{X) = g(X) - p (A) + q(A), 
and vanishes at the point A = u c . It can be expressed explicitly as 



(3.17) 



(3.18) 



4>(X; x, t) = -y/u c - XF(u c ; x, t) + -(u c - X) 2 F'(u c ) x, t) 

- ±( Uc -X)l F "(u c ;x,t) -±£F"'(ti;x,t)(Z-X)ldZ, (3.19) 

where F(A; x, t) = -x + 6At + f L (X). 

The sign of the real part of eft is such that the exponentials in the jump matrices 
decay uniformly as e — > on the jump contour, except in the vicinity of u c , where <f> 
behaves like 

0(A;x c ,t c ) = -c(u c - A) 7/2 (l + C(A-n c )), as A -> u c . (3.20) 
The constant c is strictly positive whenever the generic condition f'l'(u c ) 7^ holds. 
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Global parametrix 



Ignoring a small neighborhood of u c and the corrections to the leading order behavior 
as e — > 0, one arrives at a RH problem with jump only on (u c , +00): the jump on (u c , 0) 

'0 1 N 



is i(T\ and the jump on (0, +00) is o\, where o\ 
solved explicitly, indeed the function 

PH(A) = (_A)V4( Uc _ A )-<r 3 /4 P_ 1 



1 



. This RH problem can be 



(3.21) 



satisfies the required jump relations and as A 
tion (c) in the RH problem for S. 



00, it behaves as prescribed by condi- 



Local parametrix 

Since we have neglected a neighborhood of u c so far, where the jump matrix vs is not 
uniformly close to the identity matrix, it is necessary to construct a local parametrix in 
the vicinity of u c , which satisfies approximately the same jump conditions than S, and 
which 'matches' in an appropriate way with the global parametrix. This construction 
is rather technical but the parametrix can be constructed explicitly in terms of 1 J/, the 
solution to the Pf RH problem discussed in Section 12.11 It appears that this is true 
independently of the precise form of the function (p. Only the local behavior of 4> near 
u c matters: whenever ^(A) ~ —c(u c — A) 7 / 2 with c > as A — > u c , the local parametrix 
can be constructed in terms of *S?. The power of vanishing 7/2 can be identified with the 
power 7/2 in the asymptotic behavior for in (|2.4p and (|2.5p . This is a first indication 
that the Pf solution U (X, T) plays a role in much more general settings. Indeed if we 
would consider the RH problem for M, but with more complicated time-dependence of 
the reflection coefficient (for example corresponding to another equation than KdV), 
the local parametrix would still be built out of ^f. This supports, from the RH point 
of view, the conjecture that U (X, T) is in some sense universal. 

Remark 3.1 If one wants to perform an asymptotic analysis of the RH problem for 
M away from the critical point (x c ,t c ) as in [12] (for example for t < t c ), the main 
difference is that 0(A) ~ —c(u — A) 3 / 2 , so <p vanishes at a slower rate. The local 
parametrix can in this case be constructed in terms of the Airy function. 
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Asymptotics for M 

Multiplying S by the inverse of the local parametrix near u c and by the inverse of the 
global parametrix elsewhere, one obtains a function R which solves a so-called small- 
norm RH problem: one can show that R tends to the identity matrix in the double 
scaling limit where e — > 0, x — > x c , t — > t c at appropriate speeds. One can even obtain 
the next terms in the asymptotic series for R, which are of order e 1//7 and order e 2 / 7 , and 
which can be expressed explicitly in terms of U(X,T). Inverting the transformations 
S h-> R, T h-> S, and M h- > T, this produces an asymptotic expansion for M and, by 
(13. 8p . for the KdV solution u. In particular the method described here allows one to 
prove that the KdV solution has an expansion of the form (jl.5p . 



Remark 3.2 The above RH analysis works fine for initial data that are such that the 
initial reflection coefficient is analytic and decays sufficiently fast at infinity. If the initial 
reflection coefficient would be C 1 or C 2 instead of analytic, the usual RH techniques 
cannot be applied, but one may hope that recently developed dbar techniques [15j SI] 
could help to get around this problem. 

4 The Riemann-Hilbert approach for other equations 

4.1 The KdV hierarchy 

The equations in the KdV hierarchy can be written as a Lax equation 

L t = [L,A m ], (4.1) 

where L is (as for the KdV equation) the Schrodinger operator L = e 2 9 2 +u and where 
A m is an antisymmetric higher order operator with leading order (_1 ) m + 1 4 m e 2m Q2m+i^ 
The lower order terms of A m are determined by the requirement that [L, A m ] is an 
operator of multiplication with a function depending on u, u x , . . . , U(2 m +\)x- For m = 1, 
the operator A was given in (13.1(1 . As explained in Section 13.1} the Jost solutions to 
the Schrodinger equation can be used to construct a RH problem. The only difference 
here is the time evolution of the reflection coefficient, which is now given by 

2 4 m ? 2m+l 

r(X;t,e) = r (X;e)e— (4.2) 

With this modification, the RH problem for M given in Section 13.11 still characterizes 
the solution to the Cauchy problem for the m-th equation in the KdV hierarchy which 
is given by fl3.8f) . In the asymptotic analysis of the RH problem, the (/-function has to 
be modified and the function (j) will be somewhat different, but the main issue is that 
it will still behave like in (13.2Q[) locally near u c . Details can be found in [7]. 

4.2 The Camassa-Holm equation 

The Camassa-Holm equation (|1.14p is the compatibility condition of the Lax equations 
« 2 /rr = ( \ - A(u - e 2 ti„ + k)J / (4.3) 
/t __a + „) + (4.4) 
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This fact has been used in [J] to construct a RH problem that characterizes solutions to 
the Camassa-Holm equation. The RH solution is, in a somewhat similar manner as for 
the KdV problem, constructed using fundamental solutions to the spatial Lax equation 
(|4.3p . and the jump matrix contains a reflection coefficient for this equation. Equation 
(|4.3p is related to the Schrodinger equation through the unitary Liouville transform |3j. 
We present the RH problem in the variable z, where z 2 = —A. 

RH problem for M CH 

(a) M CH : C \ R -> C lx2 is meromorphic for z G C\R, 

(b) M has continuous boundary conditions M±(\) for z£t that satisfy the jump 
conditions 

Mr(-)=^(-)(i^g;C r ^f~ 2VyZ ) for ZEE, (4.5) 

(c) as z — > oo, we have 

M CH (z) = (l l)+o(l). (4.6) 

Again, in general there is a finite number of poles where one has to add residue con- 
ditions to ensure the uniqueness of a RH solution. In the simplest case there are no 
eigenvalues and M = M (z; y, t, e) is holomorphic in C \ M. 
The CH solution u(x,t,e) is recovered from the relation 



u(M ' f)= lV ( 2 ;!/ ' (,£ ' 



where y(x, t) is given by 



(4.7) 



y=y(x,t) 



The time evolution for the reflection coefficient is explicitly given: 

f Aitz 1 \ 

r(z; t, e) = r {z; e) exp ( — - - J . (4.8) 

This RH problem has been used to obtain large time asymptotics for CH solutions 
[31 d] . In the small dispersion setting, there is an additional complication because the 
reflection coefficient depends in a nontrivial way on the small parameter e. If one makes 
the ansatz that the initial reflection coefficient is analytic in z and behaves like 

/ 2i \ 

ro(z; e) ~ i exp ( p(z) J , as e — > (4-9) 

for a suitable function p, one could proceed in a similar manner as for the KdV RH 
problem, with the construction of a (7-function and the opening of the lens. For an 
appropriate choice of n, the ^-function should satisfy the conditions 

(a) g is analytic in C \ [— u, u], 
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(b) g has the jump relations 



g+(z) +g^(z) -2p(z) + 2f3(z;y,t) = 0, for A G (-u,u), 
(c) g(z) — > as z — > oo. 



Here /3 is given by /3(z; y, t) = —z yy — l+iz z tj • After the construction of the g-function, 

the lens can be opened as in the KdV case. The jump matrices are then supposed to 
tend to the identity matrix except near ±u. One expects that the analogue of 4> will 
behave like c\u =F z\ 7 ^ 2 near ±u for x = x c , t = t c , and that the local parametrix can 
be constructed in terms of the RH problem for but many technical issues have to 
be solved in order to proceed with this method. Another nontrivial problem one has 
to overcome to prove rigorously the universality conjecture for Camassa-Holm, is to 
identify a class of initial data such that r$ is analytic and such that (|4.9p holds. 



4.3 The de-focusing NLS equation 

The NLS equation (|1.16p is the compatibility condition of the Lax pair equations 

* x = A$, 9 t = BV, (4.10) 

where 

A = -(~1 Z l), (4.11) 



e \ u iz r 

B= \_(-2iz^-iH ?*u + «0 (4.12) 
e \ 2zu — iu x 2iz A + i\u\ .1 

The spatial Lax equation is the Zakharov-Shabat equation, and a scattering theory 
for it has been developed in [H HU EI]- Using special solutions to it, one constructs a 
solution to the following matrix RH problem. 

RH problem for M NLS 

(a) M NLS : C \ M -> C 2x2 is meromorphic for z G C\R, 

(b) M NLS has continuous boundary conditions M± LS (z) for z£M that satisfy the 
jump conditions 

M? LS (z) = M NLS (z)( 1 -^ Z l\f 2 -f{z;t,e)e-^\ for z G M, (4.13) 
y r(z; t, e)e « IZ 1 J 

(c) as z — > oo, we have 

M NLS (z;x,t,e) = I + 0(z- 1 ). (4.14) 

The NLS solution u(x, t, e) is recovered by 

u(x, t, e) = -2 lim zM^ 2 LS (z; x, t, e) (4.15) 

The reflection coefficient r(z;t,e) has a time evolution given by 

r(z;t,e) = r (z;e)e^ tz2 . (4.16) 
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If one makes the ansatz ro(z; e) ~ e~~^ ^ for a suitable analytic function /o, one can 
start performing a Deift/Zhou steepest descent analysis for this RH problem that should 
be somewhat similar to the analysis for the KdV equation. An important question is 
then to find a class of e-independent initial data for the de-focusing NLS equation that 
correspond to such reflection coefficients. 



4.4 The focusing NLS equation 

The RH problem for the focusing NLS equation is the same as for de- focusing NLS, 
but with the modified jump relation 

M? LS (z)=M» LS (z)( 1 ^ r{ f^ r(z;t,e)e-^*\ ^ g £ R _ ^ 
y r[z y t, e)e * 1 J 

Although this may seem to be a minor change compared to the de-focusing case, it has 
important consequences. Small dispersion asymptotics for the focusing NLS equation 
away from critical points have been studied for example in [361 1481 14"9]. In [2], critical 
behavior was studied under the assumption that the initial reflection coefficient is of 
the form ro(z) = e~~^°^ z \ where /o needs to satisfy some assumptions. Because of 
the sign changes in the jump matrix, one cannot construct a suitable (/-function which 
is analytic except on an interval of the real line, instead g will have its branch cut 
on a curve in the complex plane. Near two special points ati,ct2, the analogue of the 
function <f> we described in the KdV case behaves like <p(z) ~ c(z — ctj) 5 / 2 . This power 
of vanishing can be identified with the power 5/2 in (12. lip , and a local parametrix 
can be constructed in terms of the RH solution <3?. At least this is true for values of 
x, t, e that are mapped to a complex number Z such that the RH problem for $ with 
parameter Z is solvable (or in other words values of Z where the tritronquee solution 
Q has no pole). The results in [2] confirm the conjectured behavior from [21 j . but two 
questions remain unanswered. First the absence of poles for the tritronquee solution 
Q(Z) to Painleve I in the sector | arg Z\ < ^ has not been proved yet. Secondly, given 
initial data for the focusing NLS equation, semiclassical asymptotics for the reflection 
coefficient (which corresponds to a non-self-adjoint Zakharov-Shabat operator) have 
not been obtained. In [2], an ansatz has been made about the small e-behavior of the 
reflection coefficient, and an asymptotic analysis of the RH problem has been done for 
such reflection coefficients. This implies a proof of the conjectured behavior for initial 
data that correspond to a reflection coefficient of the proposed form, but it is not known 
that the reflection coefficient corresponding to appropriate e-independent initial data 
is of the prescribed form. 



5 Conclusion and outlook 

For Hamiltonian perturbations of hyperbolic equations or systems which can be 'solved' 
by inverse scattering, or in other words which can be characterized in terms of a RH 
problem, one can hope at this point to prove that they obey Dubrovins universality 
conjecture. The Deift/Zhou steepest descent method allows one to analyze the RH 
problem asymptotically in the small dispersion limit e — > 0, although several of the 
RH transformations are considerably different for different equations. In particular the 
g-function mechanism has to be handled separately for each equation. This means 
that it is not obvious at all how one can prove the universality simultaneously for a 
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whole class of equations. Another issue is that one needs to have detailed asymptotic 
information about the reflection coefficient in the semiclassical limit e — > before one 
can start with the asymptotic analysis of the RH problem. General equations of the 
form (jl.3p cannot be solved using inverse scattering; for such 'non-integrable' equations 
there are no rigorous results about the asymptotic behavior available, even existence 
of solutions has not been established in general. 

In the elliptic case, much less is clearly understood as we explained in the previ- 
ous section. Although the tritronquee solution to Painleve I has a long history, the 
conjecture about its poles remains open. Another challenge is to obtain semiclassi- 
cal asymptotics for the reflection coefficient related to the Zakharov-Shabat equation. 
Apart from the focusing NLS equation, no rigorous results have been obtained concern- 
ing critical behavior in the elliptic case. 
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